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Path integrals have played a fundamental role in emphasizing the profound 

■ analogies between Quantum Field Theory (QFT), and Classical as well as Quan- 
Q\ \ turn Statistical Physics. Ideas coming from Statistical Physics have then led to a 

deeper understanding of Quantum Field Theory and open the way for a wealth 
of non-perturbative methods. Conversely QFT methods are become essential 

■ for the description of the phase transitions and critical phenomena beyond mean 
[ field theory. This is the point we want to illustrate here. We therefore review 

""p^ ■ the methods, based on renormalized </>t quantum field theory and renormaliza- 

tion group, which have led to an accurate determination of critical exponents of 
q^. the iV-vector model, and more recently of the equation of state of the 3D Ising 

(D [ model. The starting point is the perturbative expansion for RG functions or the 

,^Zh • effective potential to the order presently available. Perturbation theory is known 

to be divergent and its divergence has been related to instanton contributions. 

• »— j ■ This has allowed to characterize the large order behaviour of perturbation series, 

an information that can be used to efficiently "sum" them. Practical summa- 
tion methods based on Borel transformation and conformal mapping have been 
developed, leading to the most accurate results available probing field theory in 
a non-perturbative regime. We illustrate the methods with a short discussion of 
the scaling equation of state of the 3D Ising model [[jj. Compared to exponents 
its determination involves a few additional (non-trivial) technical steps, like the 
use of the parametric representation, and the order dependent mapping method. 
A general reference on the subject is 

J. Zinn-Justin, 1989, Quantum Field Theory and Critical Phenomena, in par- 
ticular chap. 28 of third ed., Clarendon Press (Oxford 1989, third ed. 1996). 
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1 Introduction 

Second order phase transitions are continuous phase transitions where the cor- 
relation length diverges. Renormalization group (RG) arguments || , as well as 
an analysis, near dimension four, of the most divergent terms appearing in the 
expansion around mean field theory ||, indicate that such transitions present 
universal features, i.e. features independent to a large extend from the details of 
the microscopic dynamics. Moreover all universal quantities can be calculated 
from renormalizable or super-renormalizable quantum field theories. For an 
important class of physical systems and models (with short range interactions) 
one is led to a </> 4 -like euclidean field theory with O(N) symmetry. Among 
those let us mention statistical properties of polymers, liquid-vapour and bi- 
nary mixtures transitions, superfluid Helium, ferromagnets... We explain here 
how critical exponents and other universal quantities have been calculated with 
field theory techniques. To simplify notation we concentrate on the universality 
class of the Ising model (models with Z2 symmetry). 

The effective quantum field theory. The relevant field theory action TC ((f)) is 



where r plays the role of the temperature. The mass parameter A corresponds 
to the inverse microscopic scale and also appears as a cut-off in the Feynman 
diagrams of the perturbative expansion. 

For some value r c the correlation length diverges (the physical mass van- 
ishes), and the critical domain corresponds to \t = r — r c \ < A 2 . The study of 
the critical domain reduces to the study of the large cut-off behaviour, i.e. to 
renormalization theory and the corresponding renormalization group. However 
one essential feature of the action distinguishes it from field theory in the form 
it was traditionally presented in particle physics: the dependence on A of the 
coefficient of (p 4 is given a priori. In particular in the dimensions of interest, 
d < 4, the "bare" coupling constant diverges, though the field theory, being 
super-renormalizable, requires only a mass renormalization. 

To circumvent the problem of the large coupling constant the famous Wilson- 
Fisher ^-expansion Q has been invented. The dimension d is considered as a 
continuous variable. Setting d = 4 — e one expands both in go and e. Diver- 
gences then behave like in four dimensions, they are only logarithmic and can 
be dealt with. Moreover it is possible to study directly the massless (or critical) 
theory. 

Later it has been proposed by Parisi ||, to work at fixed dimension d < 4, 
in the massive theory (the massless theory is IR divergent). One motivation 





2 



for trying such an approach is of practical nature: it is easier to calculate 
Feynman diagrams in dimension three than in generic dimensions, and thus 
more perturbative orders are available. 

The large A limit then is taken first at uq = groA 4_d fixed. This implies 
that one first tunes the initial parameters of the model to remain artificially 
close to the unstable uq = gaussian fixed point: When the correlation length £ 
increases near T c one decreases the dimensionless bare quantity go as fi^ 4 ^ oc 
1/£A. Finally one takes the infinite uq limit. One then is confronted with a 
serious technical problem: perturbation theory is finite but one is interested in 
the infinite coupling limit. 

One thus introduces a field renormalization, (f> = Z 1 / 2 ^, and a renormalized 
dimensionless coupling constant g as in four dimensions. They are implicitly 
defined by the renormalization conditions for the 0r 1PI correlation functions: 

rg (p; m,g) = m 2 + p 2 + O (p 4 ) , ( Pl = 0; m, g) = m^ d g . (1.2) 

The role of renormalizations, however, is here different. When the corre- 
sponding Callan-Symanzik /3-function has an IR stable zero, /3(g*) = with 
u = P'(g*) > 0, then the new coupling g has a finite limit g = g* when the 
initial coupling constant uq becomes large. Thus the renormalized coupling g 
is a more suitable expansion parameter than uq. 

Note that the mass parameter m, which is proportional to the physical mass, 
or inverse correlation length, of the high temperature phase, behaves for t = 
r — r c — > 0-|- as m oc f , where v is the correlation length exponent. 

In contrast with the ^-expansion however, at fixed dimension three or two 
one has no small parameter. Therefore accurate determinations of g* and all 
other physical quantities depends on the analytic properties of the series, in 
addition to the number of terms available. A semi-classical analysis, based on 
instanton calculus, unfortunately indicates that perturbation theory in </> 4 field 
theory is divergent. Therefore to extract any information from perturbation 
theory a summation method is required. 

2 Critical exponents. Borel transformation and mapping 

The most studied quantities are critical exponents, because they are easier to 
calculate. They have been extensively used to compare RG predictions with 
other results (experiments, high or low temperature series expansion, Monte- 
Carlo simulations). The first accurate values of the exponents of the O(N) 
symmetric N- vector model have been reported in ref. using six- loop series 
for RG functions [0 . Perturbative series have been summed using Borel trans- 
formation and conformal mapping. The same ideas have later been applied to 
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the £-expansion when five loop series have become available, and recently to 
the equation of state. With time the method has been refined and the efficiency 
improved by various tricks but the basic principles have not changed. 

Borel transformation and conformal mapping. Let R(g) be a quantity given 
by a perturbation series 

R{g) = Y,Rk9 k . (2.1) 

Large order behaviour analysis (instantons) || teaches us that, in the </> 4 field 
theory, Rk behaves like k s (—a) k kl for k large. The value of a > has been 
determined numerically. One thus introduces B(g), the Borel transform of 
R(g), which is defined by 

B(g) = y £(Rk/kl)g k . (2.2) 

The function B(g) is analytic at least in a circle with the singularity closest to 
the origin located at z = — 1/a. Unlike R(g), B(g) is determined by its series 
expansion. In the sense of formal series, R(g) can be recovered from 

POO 

R(g)= / e-*fl(0t)d*, (2.3) 
Jo 

However, for relation ( |2.3|) to make sense as a relation between functions, and 
not only between formal series, one must know B(g) on the whole real positive 
axis. This implies that B(g) must be analytic near the axis, a result proven 
rigorously |J. Moreover it is necessary to continue analytically the function 
from the circle to the real positive axis. Consideration of general instanton 
contributions suggests that B(g) actually is analytic in the cut- plane. Therefore 
the analytic continuation can be obtained from a conformal map of the cut- 
plane onto a circle: 

z u(z) = az/(y/l + za + l) 2 . (2.4) 
The function R(g) is then given by the new, hopefully convergent, expansion 

POO 

R{g) = Y^B k e"* [«((/*)]* d*. (2.5) 
Jo 

Exponents. The values of critical exponents obtained from field theory have 
remained after about twenty years among the most accurate determinations. 
Only recently have consistent, but significantly more accurate, experimental 
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Critical exponents of the O(N) models from d = 3 expansion ]?!}/ . 



N 





1 


2 


3 


* 

9 


26.63 ±0.11 


23.64 ±0.07 


21.16 ±0.05 


19.06 ±0.05 


7 


1.1596 ±0.0020 


1.2396 ± 0.0013 


1.3169 ±0.0020 


1.3895 ± 0.0050 


v 


0.5882 ±0.0011 


0.6304 ± 0.0013 


0.6703 ±0.0015 


0.7073 ±0.0035 


V 


0.0284 ± 0.0025 


0.0335 ± 0.0025 


0.0354 ±0.0025 


0.0355 ± 0.0025 


(3 


0.3024 ± 0.0008 


0.3258 ± 0.0014 


0.3470 ±0.0016 


0.3662 ± 0.0025 


u> 


0.812 ±0.016 


0.799 ±0.011 


0.789 ±0.011 


0.782 ±0.0013 



results been reported in low gravity superfluid experiments ||10|| . Also various 
numerical simulations jll] and high temperature expansions on the lattice have 
claimed similar accuracies. 

The values of critical exponents have recently been updated |T^] because 
seven-loop terms have been obtained for two of the three RG functions. Some 
results are displayed above. The main improvements concern the exponent rj 
which was poorly determined, and the lower value of 7 for N = (polymers). 



3 3D Ising model: the scaling equation of state 



Let us first recall a few properties of the equation of the state in the critical 
domain, in the specific case N = 1 (Ising-like systems), at d = 3. 

The equation of state is the relation between magnetic field H , magnetization 
M = {(f)) (the "bare" field expectation value) and the temperature which is 
represented by the parameter t = r — r c oc T — T c . It is related to the free energy 
per unit volume, in field theory language the generating functional T(</>) of 1PI 
correlation functions restricted to constant fields, i.e the effective potential V, 
V(M) = r(M)/vol., by H = dV/dM. In the critical domain the equation of 
state has Widom's scaling form 



H(M, t) = M^it/M 1 / 13 ), (3.1) 

a form initially conjectured and which renormalization group has justified. 

One property of the function H(M, t) which plays an essential role in the 
analysis is Griffith's analyticity: it is regular at t = for M > fixed, and 
simultaneously it is regular at M = for t > fixed. 

In the framework of the ^-expansion, the function f(x) has been determined 
up to order e 2 for the general O(N) model, |TB|] and order e 3 for N = 1 |JL4|| . 
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The calculations presented here are performed within the framework of the 0| 
massive field theory renormalized at zero momentum (eq. ( |1 . 2| ) ) . Five loop se- 
ries for the effective potential have been reported [|15| . The conditions (|1.2| ) im- 
ply that the effective potential V expressed in terms of the expectation value 
of the renormalized field ip = (</>r), has a small <p expansion of the form 

V(<p) - V(0) = ±mV + h m 9^ + O (</) = (m s /g)V(z, g), (3.2) 

where z is a dimensionless variable z = ip^/g/m. For g = g* 

^«M/Vm^ocM/m (1+,I) / 2 «M^ ,3 . (3.3) 

The equation of state is related to the derivative F of the reduced effective 
potential V with respect to z 

H<xt? s F(z), F(z)=F(z, 9 n = dVi *f ) . (3.4) 

The problem of the low temperature phase. To determine the equation of 
state in the whole physical range a new problem arises. In this framework it is 
difficult to calculate physical quantities in the ordered phase because the theory 
is parametrized in terms of the disordered phase correlation length £ = m _1 oc 
t~ v which is singular at T c (as well as all correlation functions normalized as 
in ( |1.2[ )). In the limit m — > 0, at ip fixed, z — > oo as seen in eq. ( |3.3[ ). In this 
limit from eq. ( |3.1|) one finds 

H(M, t = 0) oc M s => F(z) oc z s . (3.5) 

The perturbative expansion of the scaling equation of state leads to an expres- 
sion only adequate for the description of the disordered phase. 

In the case of the £-expansion the scaling relations (and thus the limiting 
behaviour (|3.5| )) are exactly satisfied order by order. Moreover the change 
to the variable x oc z~ x l^ (more appropriate for the regime t — ^ 0) gives an 
expression for f(x) oc F(x~ /3 )x^ 3S that is explicitly regular in x = (Griffith's 
analyticity). Still, even there a numerical problem arises when e = 1 is set. 

In the case of fixed dimension perturbation theory, instead, because IR scal- 
ing is obtained only for g = g* and not for generic values of g, scaling properties 
are not satisfied at any finite order in g. 

Several approaches can be used to deal with the problem of continuation 
to the ordered phase. A rather powerful method, motivated by the results 
obtained within the e-expansion scheme, is based on the parametric represen- 
tation. 
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4 Parametric representation of the equation of state and ODM 

Both the scaling and regularity properties of the equation of state can be more 
easily expressed by parametrizing it in terms of two new variables R and 6, 
setting: 

M = m R p 9, t = R(l-9 2 ) 1 H = h R l3S h(e) , (4.1) 

where ho : mo are normalization constants. Then the function h{6) is an odd 
function of 6 which from Griffith's analyticity is regular near 6 = 1, which is z 
large, and near 6 = which is z small. It vanishes for 6 = 6q which corresponds 
to the coexistence curve H = 0, T < T c . In terms of the scaling variable z used 
previously one finds 

z = pd/(l-6 2 f, h(6) = (l-6Y s F(z(6)), (4.2) 

where p is an arbitrary parameter. Note if we know a few terms of the expansion 
of F(z) in powers of z we know the same number of terms in the expansion 
of h(6) in powers of 6. But because h is a more regular function, the latter 
expansion has a much larger domain of validity. 

From the parametric representation of the equation of state it is then possible 
to derive a representation for the singular part of the free energy per unit 
volume as well as various universal ratios of amplitudes. 

Order dependent mapping (ODM). In the framework outlined before, the 
approximate h(6) that one obtains by summing perturbation theory at fixed 
dimension, is still not regular. The terms singular at 6 = 1, generated by the 
mapping (|4.2|), do not cancel exactly due to summation errors. The last step 
thus is to Taylor expand the approximate expression for h(6) around 6 = and 
to truncate the expansion, enforcing in this way regularity. A question then 
arises, to which order in 6 should one expand? Since the coefficients of the 6 
expansion are in one to one correspondence with the coefficients of the small z 
expansion of the function F(z), the maximal power of 6 in h(8), should be equal 
to the maximal power of z whose coefficient can be determined with reasonable 
accuracy. As noted before, although the small z expansion of F(z) at each 
finite loop order in g contains an infinite number of terms, the evaluation of 
the coefficients of the higher powers of z is increasingly difficult. 

Therefore one has to ensure the fastest possible convergence of the small 6 
expansion. For this purpose one uses the freedom in the choice of the parameter 
p in eq. Q4.2| ): one determines p by minimizing the last term in the truncated 
small 6 expansion, thus increasing the importance of small powers of 6 which 
are more accurately calculated. This is nothing but the application to this 
particular example of the series summation method based on ODM [ 117) . 
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5 Numerical results 

One first determines the first coefficients -F2Z+1 of the small z (small field) 
expansion of the function F(z) as accurately as possible, using the same method 
as for exponents, i.e. Borel-Leroy transformation and conformal mapping. One 
finds 

F 5 = 0.01711 ±0.00007, F 7 x 10 4 = 4.9 ± 0.5 , F 9 x 10 5 = -7 ± 5 . 

One then determines by the ODM method the parameter p and the function 
h{9) of the parametric representation, as explained before. One obtains suc- 
cessive approximations in the form of polynomials of increasing degree. At 
leading order one obtains a polynomial of degree five. It is not possible to go 
beyond hg(p) because already F 9 is poorly determined. Note that one has here 
a simple test of the relevance of the ODM method. Indeed, once h(9) is deter- 
mined, assuming the values of the critical exponents 7 and (3 one can recover 
a function F(z) which has an expansion to all orders in z. As a result one 
obtains a prediction for the coefficients F21+1 which have not yet been taken 
into account to determine h(6). The relative difference between the predicted 
values and the ones directly calculated gives an idea about the accuracy of the 
ODM method. The simplest representation of the equation of state, consistent 
with all data, is given by 

h(6) = 9- 0.76201(36) 9 s + 8.04(11) x 10" 3 9 5 , (5.1) 

(errors on the last digits in parentheses) that is obtained from p 2 = 2.8667 . 
This expression of h(9) has a zero at #0 = 1-154, which corresponds to the 
coexistence curve. The coefficient of 9 7 in eq. ( |5.1| ) is smaller than 10 -3 . Note 
that for the largest value of 9 2 which corresponds to 9 2 , the 9 5 term is still a 
small correction. 

Concluding remarks. Within the framework of renormalized quantum field 
theory and renormalization group, the presently available series allow, after 
proper summation, to determine accurately critical exponents for the A-vector 
model and the complete scaling equation of state for 3D Ising-like (A = 1) 
systems. In the latter example additional technical tools, beyond Borel sum- 
mation methods, are required in which the parametric representation plays a 
central role. From the equation of state new estimates of some amplitude ratios 
have been deduced which seem reasonably consistent with all other available 
data. 

Acknowledgments. Hospitality at the Center for Theoretical Physics, 
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